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Abstract 

We study the Q 2 dependence of the chiral-odd twist-3 distribution e(x, Q 2 ). The anoma- 
lous dimension matrix for the corresponding twist-3 operators is calculated in the one-loop 
level. This study completes the calculation of the anomalous dimension matrices for all 
the twist-3 distributions together with the known results for the other twist-3 distributions 
g2{x, Q 2 ) and Hl(x, Q 2 ). We also have confirmed that in the large N c limit the Q 2 -evolution 
of e(x, Q 2 ) is wholely governed by the lowest eigenvalue of the anomalous dimension matrix 
which takes a very simple analytic form as in the case of and Hl- 
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1 Introduction 



Ongoing plans of high energy collider experiments are going on their way toward more 
precision measurements and more variety of spin-dependent observables in many (semi) 
inclusive processes. Correspondingly, increasing attention is paid on the power corrections 
due to the higher twist-effects which represent parton correlations in the target. Recent 
measurement of the nucleon's #2 structure function by the E143 collaboration JIJ anticipates 
a forthcoming significant progress of the twist-3 physics. Under this circumstance, QCD 
study on the scale dependence of various distribution and fragmentation functions is of great 
interest . 

The nucleon has three twist-3 distributions g2, and e 0. #2 is chiral-even and the other 
two are chiral-odd. e is spin- independent and the other two are spin-dependent. Compared 
to e(x,Q 2 ), g2 and have more chance to be measured experimentally since they become 
leading contribution to proper asymmetries in the polarized deep inelastic scattering (DIS) 
and Drell-Yan processes, respectively. Their Q 2 evolution has been studied in |], |5| for #2 
and in || [7j for h^. 

Similarly to the distribution functions, there are three twist-3 fragmentation functions 
which describes hadronization processes of partons in semi-inclusive processes ||, §2(2, Q 2 ), 
iiL(z,Q 2 ) and e(z,Q 2 ). (Their naming is parallel to the corresponding distribution func- 
tions.) In the inclusive pion production in the transversely polarized DIS, the chiral-odd 
fragmentation function e of the pion appears as a leading contribution together with hi 
(twist-2) of the nucleon. Although the Q 2 evolution of the twist-2 fragmentation functions 
is known to be obtained from that of the corresponding distributions in the one-loop level 
(Gribov-Lipatov reciprocity no such relation is known for the higher twist fragmentation 
functions. 

In this paper we investigate the Q 2 evolution of e(x, Q 2 ). Theoretically, this completes 
the calculation of the anomalous dimension matrices of all the twist-3 distributions. Phe- 
nomenologically, we expect it will shed light on the Q 2 evolution of e(z, Q 2 ), anticipating the 
day when their relation is clarified. 

The outline of this paper is the following: In section 2, we briefly recall the twist-3 
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operators for e(x, Q 2 ) 0. Similarity to the hi case is reminded. In section 3, we present the 
calculation of the one-loop anomalous dimension matrix for e(x,Q 2 ). The method follows 
that of 0. In section 4, the Q 2 evolution of e(x, Q 2 ) in the large N c limit is discussed. This 
part is a recapitulation of J7| in our language. Section 5 is devoted to a brief summary of 
our results. Appendix contains the contributions from each one-loop Feynman diagram. 

2 Twist-3 operators for e(x, Q 2 ) 

The chiral-odd twist-3 distribution function e(x, Q 2 ) is defined by the relation g] 

/ ^ e ^(P\m^n)\ Q \P) = 2Me(x,Q 2 ) } (2.1) 

J 271 

where \P) is the nucleon (mass M) state with momentum P. Two light-like vectors p and n 
defined by the relation, P = p+ ^f-n, p 2 = n 2 = 0, p-n = 1, specify the Lorentz frame of the 
system. Gauge-link operators are implicit in (|2.1| ). Taking the moments of ( |2.1|) with respect 



to x, one can express the moments of e(x, Q 2 ) in terms of the nucleon matrix elements of 
the twist-3 operators V^ 1 

M n [e(Q 2 )} = e n (Q 2 ), (2.2) 
(PS\V^ 2 ~^ n \PS) = 2e n M(P fll P f * 2 ...p^- traces), (2.3) 
ynw-nn = SnTpiD^iD 1 * 2 ■■■iD^ip- traces, (2.4) 

where M. n [e(Q 2 )] = j\ dxx n e(x, Q 2 ) and the covariant derivative = 3^ — igA^ restores 
the gauge invariance. As in (|2.3|) and (|2.4j ), we often suppress the explicit scale dependence. 
Following a common wisdom we introduce a null vector A^ (A 2 = 0) to kill the trace 
terms of V^ 1 '"^" and introduce the notation V n ■ A = p 1 '"'* B A w ■ ■ ■ A^ for convenience, 
and similarly for other operators in the following. Using the relation = | {p,7 M } and 
[D^ D v \ = —igG^u with the gluon field strength G^ u , V n ■ A can be recast into the following 
form: 

n 

V n ■ A = U n ,r A + N n ■ A + E n ■ A, (2.5) 

1=2 

where 

Up'""" = -Sn^a^HD^-'-gG^-'-iD^i}- traces, (2.6) 



N n~t*» = Snrriqijj^iD^ ■ ■■iD' Xn ijj -traces, (2.7) 
1 
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^(ip - m q )^HD^ ■ ■ ■ iD^ip + tyfHD** ■ ■ ■ iD^ n (ip - m q )ip 



traces. 
(2.8) 



Ui contains G^ v explicitly, which indicates that e(x) represents the quark-gluon correlations 
in the nucleon. N is the quark mass times the twist- 2 operator which contributes to the 
moments of fi structure function familiar in the spin- averaged DIS. E is the EOM (equation 
of motion) operator which vanishes by use of the QCD equation-of-motion. Although the 
physical matrix elements of EOM operators vanish, one needs to take into account the mixing 
with E to carry out the renormalization of U and N , as is discussed in |5], |6| in the context 



of the renormalization of g% and Hl- From ( |2.5| ) one sees that Ui appears in the form of 

+ l) . (2.9) 



pA t l'"A t n TTt J -l---l J -n , T r/il • • -fin I 1 9 



By this combination, R n i has a definite charge conjugation (— l) n . Readers may recall the 
similarity between the present {Ui, R n j} and {61, R n j} which appeared in h L P, [|. In fact 
presence of 75 in 61 is the mere difference from U\. R n j of hi is defined as n _i +2 — Q\ and 
has a charge conjugation (— l) n+1 which is opposite to the above R n j in Q2.9p. We will see 
the similarity in the renormalization constants between and e in the following. 

3 Q 2 evolution of e(x, Q 2 ) 

For the renormalization of e(x), we closely follow the method of |§ which discussed the 
renormalization of Hl- So we omit the detail in the following. Owing to the chiral-odd 
nature, e(x) does not mix with the gluon-distribution. We choose Rn,i (1 = 2,..., | +1), 
E, N as a basis of the operators. For the renormalization of Rnj, we calculate the one- loop 
correction to the three-point function F^p, q,p — q) defined by 

F^p, q, fc)(27r) 4 <5 4 (p + k - q)G(p)G(q)D(k) 

d 4 xd 4 yd 4 ze ipx e- lli ye lkz (T{0^(x)^(y)A,(z)}), (O = R nh N, E) (3.1) 



where G and D are, respectively, the quark and gluon propagators. (We suppressed the 
Lorentz and spinor indices for simplicity.) In order to take into account the mixing with the 



EOM operator properly, we use off-shell kinematics for the external lines. The calculation 
is done with the Feynman gauge for the gluon propagator (of course the result should be 
independent of the gauge) and the MS scheme is adopted with the dimensional regularization. 
The three-point basic vertices shown in Fig. 1(a) for R n j ■ A, E n • A, N n ■ A are calculated 
to be 



-i^a aX A x (p n - l q 1 - 2 +tf- 2 q n - l )(-kg alt + k a AJt a , 



(3.2) 
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3=1 



m. 



3=2 
n 

3=2 



t a . 



(3.3) 
(3.4) 



where k = q — p, p = p- A for an arbitrary four vector p and t a is the color matrix normalized 
as Tr(t a t b ) = \$ ab - One-loop diagrams for the above are shown in Fig. 2. In calculating 
these diagrams, we need a Feynman rule for the four-point basic vertex of R Ut i ■ A shown in 
Fig. 1(b). It is given by 



9_ 
2 



n-lA-2 



f abc t c a^A x A,g au p n ~ l q 

n—l+l 

*A a j-a+b 



— /cr~ A A\t a t b P J - 2 A,(p + kr- l+1 ^(-k'g ai/ + k' a A u )t 

3=2 



-ta 



'A x t a t b Yl P n ~ l K(p + ky- n+l -\-kg a , + KA^ 



j=n— 1+3 

+{ji <-n/,k <-> k',a <->b)] + (I —> n — 1 + 2). 



(3.5) 



For the actual calculation of the one-loop diagrams, we introduce a vector with the 
condition Q ■ A = 0, and consider F^Q^ ||. This way mixing with the gauge noninvariant 
EOM operators can be avoided and the calculation is greatly simplified. (See || for the 
detail.) The contraction of „ an d ^nia with fi M leads to 



= l {n a a^A x (q-p)(p n - l t 2 +P l - 2 r- l )t a , 



(3.6) 
(3.7) 



The mixing coefficients between R n ^ and {N, E}, Zie and Zin, can also be obtained from the 
one-loop correction to the two point function shown in Fig. 3 (b), (c), giving a consistency 
check. The two point basic vertices for E and N are 

1 



P(2) 



-p n - 1 (A^+^A-2m,A), 



m q p 



sn-l , 



(3.8) 
(3.9) 



Using these Feynman rules for the vertices, one can calculate the one-loop diagrams 
shown in Figs. 2 and 3. Actual calculation is very tedious and will be described in the 
Appendix in detail. Taking into account the wave function renormalization for the fields and 
the renormalization for m q and g which explicitly appear in the vertices |J , we eventually 
obtained the renormalization constants Zij among R n j, N and E in the following matrix 
form: 



I Zi m 

E* = Z EE (fi) 
\N*J V Z nn (jj.) 



Zij can be expressed as 



Sa + 



-2—Y.. 

16tt% 13 



hJ 



E n {ii) 



Z, m 





"n" 












.2. 





1 .(3.10) 



+ 1,E,N 



(3.ir 



with s = -=- (d is the space-time dimension in the dimensional regularization) and the 
constants given below. To write down Y^- we define a symbol (~) as (|) = [|] = ~ for 
an even n and (|) = [|] + 1 = for an odd n. Then is given as follows: 



Y lr 



a 



G 



(m — l)(m + /) (m — l)(m + n — Z + 2) 



2(Z-m)[Z-l] 2 2(n-Z-m + 2)[n-Z + 1] 2 Z ra-Z + 2 

Z-l)(n-Z-m + 3) n-Z + 4 



(Z-3)(Z + l-m) Z + 2 (n 



2[Z - 1] 3 
+ (2C F - C G 



2[Z - 1], 



2[n - I + 1] 



2[n-Z + l] 2 



(-1) 



n—l—m 



i.C, 



n-l^m-2 



+ 



■iy- m ,_ 2 c,. 



-2(-l)" 
(2 < Z < 



(n — Z + 2 — m) n-m+lQ-l (^ — m ) n-m+lCl-m 
n-i+l'-'m-l 

t "F^]r + [n-Z+l] 3/ 
+ 1, 2 < m < Z - 1). 



(3.12) 
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Y u = C G 



o __ q 1 1 | (/-l)(n + 2) 

i_1 ™~ m 21 2{n - I + 2) 2(n - 21 + 2) [n - I + 1] 2 



n-l+2 I 

1-3 1 + 2 (n-l-l)(n-2l + 4) n-l + A 

+ 77T, 7T~ + 77T, 7T- + ~ ^ A 71 " + 



2[l - 1] 3 2[l - 1] 2 
-1)' 



2[n - Z + 1] 3 [n - I + 1] 2 



+2(C G - 2C F ) 

-(-ir 



1 ( — l) n + n _; + iC/_i 



_[Z-1] 3 [n-Z + l] 3 



Z- 1 



2(n-Z + l)(n-2Z + 2) 
-C F (25,_! + 2S n . l+1 - 3) ^2 < Z < 



(3.13) 



— 2Sr 

2 n + 2 



f-2 f+3 



+ (2C F - C G )4- 



a 



G 



+C f (3-4:Sr) (for even n), 



(2n — Z — m + 4) (n — m + 1) (m — 1) (m + n — Z + 2) 



2(m - l)[n - I + 1] 2 2(n - m - Z + 2)[n - Z + 1] 2 



2 (n-Z-l)(n-2Z + 4) n-Z + 4 



n-Z + 2 ' 2[n-Z + l] 3 [n - I + 1] 2 

2( — l) m { n _; + iC m _i + ( — 1)™ n-l+lC m -i} 



+(C G - 2C F 
_(_!)"»-' (- 



[n-l + 1] 3 

( — I)" n-jCm-2 



+ 



\(m - Z) m _iC m _/ (n - Z + 2 - ra) „_ m+ iCn 
(3< i + l< m <(|)), 



^ § +1 = Cg 



1 



n(3n - 2Z + 6) 



+ 2 



4 [n - Z + 1] 2 n - 2Z + 2 n-Z + 2 
1 f(n-Z-l)(f -Z + 2) i n-Z + 4 



[n - I + 1] £ 



+ 



[n - Z + 1]; 



+ (2Cf_Cg) ([n-Z + l] 3 



Cll + 



f-Z + 1 ^Q-x . 



(3.14) 



(3.15) 
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for even n, 2 < I < — , . 



(3.16) 



Y, 



IE 



-2C7 



2 < / < 



n 



+ 1 



(3.17) 



IN 



\l-l] 3 [n-Z + lU 



^ee — 2(1 — SnjCp, 



2 < Z < 



(3.18) 



(3.19) 



Y, 



NN 



n(n + 1) 



^3.20) 



where Cp — N £ N 1 and C*g = iV c are the Casimir operators of the gauge group SU(N C ), 
S n = Sj=i = 3 (J + 1) ■ ■ • (j + ^ — 1), and n C\ is the binomial coefficient defined as 

n Q = n\/U(n -1)1 Note that (|3~T^ ) is one half of (gl5|) with m = f + 1 for even n. With 



these Yy (i,j = 2,3,..., ^ + 1, E, N), the anomalous dimension matrix for the twist-3 
operators R n j, E and N take the form of the upper triangular matrix as 

,2 



g 



in 



-Y 

2 J W 



(3.21) 



The above results for Y can be compared with the mixing matrix X for hi given in Eqs. 
(3.14)-(3.20) of ||. Their difference comes from the opposite charge conjugation symmetry 
of R n ,i- Solving the renormalization group equation, the Q 2 evolution of R n j and iV is given 
by 



a n ,i(Q 2 



d n (Q 2 ) 



E 



'a(Q 



2\\ -Y/bo' 



m=2 

fa(Q 



a(/i 2 ) 

2\ \ ~ Y NN/bo 



I in 



'a(Q 2 

M^ 2 



-Y/bo' 



d n (fi 2 ), (3.22) 



IN 



Va(/i 2 ) J 
where a n> i and d n are defined by 

(p\K)r n ^ 2 )\p) 

(P\N^ n {fi 2 )\P) 



d n (fi 2 



(3.23) 



2a rii i(fi 2 )MS n (P f11 ...pi*- traces), 
2d n (fi 2 )MS n (P^ ...P^- traces), 



(3.24) 
(3.25) 



and b = ±±N C - \N f . 



The complete spectrum of the eigenvalues of the anomalous dimension matrix Q3.21 
(ignoring the factor g 2 /87i 2 ) is shown in Fig. 4(a) together with those for the twist-2 dis- 
tribution h\. One sees from the figure that the anomalous dimensions of e are significantly 
larger than those of h%, except for the first several moments. 

Before closing this section, we compare the Q 2 evolution among the twist-2 and -3 dis- 
tributions (at m q = 0), using the first few moments. For n = 2 and 3, only one operator 
contributes to e(x). The values of Y can be read from ( |3.13| ) and ( |3.14| ) as Y n=2 
Y n=3 = ^p. Thus the Q 2 evolution becomes 



|^ and 



M 2 \e{Q 2 )] = L & - ll ' b °M 2 k/x 2 )l , M 3 \e(Q 2 )} = L s - n / b °M 3 fe(// 



(3.26) 



where L = y 2 l . In table 1, we summarize the anomalous dimensions of the twist-2 and -3 
distributions for n — 2, 3, ignoring the common factor g 2 /8n 2 . One sees that these moments 
of e evolve slower than those of the other twist-3 distributions and are close to the twist-2 
distributions. 



n 


fi, 9i 


h 


e 


92 


h L 


2 


50 


52 


55 


77 






9 


9 


9 


9 




3 


314 


64 


73 




104 




45 


9 


9 




9 



Table 1 . The anomalous dimensions for the 2nd and 3rd moments of the nonsinglet twist-2 
and twist-3 distributions. For these moments, the twist-3 distributions receive the contribu- 
tion from only one operator, hi and g 2 denote the twist-3 parts of hi and g 2 , respectively. 
For Hl, the third moment is the lowest one, and the result is taken from 0. For g 2 , the 
result is taken from 0, and the anomalous dimensions for odd moments are not available 
because they have been discussed in the context of DIS. 



4 Large N c limit 

In it has been proved that all the (nonsinglet) twist-3 distributions g 2 , hi and e obey a 

simple GLAP equation [TO] similarly to the twist-2 distributions in the iV c — > oo limit. In this 
limit, Q 2 evolution of these distributions is completely determined by the lowest eigenvalue 
of the anomalous dimension matrix which has a simple analytic form. For e, this can be 
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checked using the results in the previous section. At N c — > oo, i.e., Cf — > N c /2, the lowest 
eigenvalue of 7 in ( |3.21 ) is given by (ignoring the factor g 2 /8ir 2 ) 

1 1 \ 



In = 2iVc S n 



4 2(n + l)/' 

As was shown in section 2, the n-th moment of e(x) can be expressed in terms of a n> i as 



(4.1) 



2 f 1 



n : even 



(4.2) 



where {• • • n : even} means this term is only for even n. We can directly check that the 
coefficients of a n i in ( |4.2|) constitute the left eigenvector of the mixing matrix Y with the 
eigenvalue —7^, i.e., 



[f]+ 1 

— — 7 n 

Z=2 



m = 2, 



+ 1 



for odd n and 



(4.3) 



L2j 1 
i=2 Z 



"7n 



m = 2, 



L2j x 



(4.4) 



for even n. This means all the right eigenvectors except the one which corresponds to 7^ are 
orthogonal to the above vector consisting of the coefficients of a ni \ in (^4.2|) . This proves that 
the Q 2 evolution of e is exactly given by 



M n \e(Q 2 )} = L^ b0 M n [e(/i 



(4.5) 



as is found in |7j by a different method. The complete spectrum of the eigenvalues of the 
anomalous dimension matrix at large iV c is given in Fig. 4(b) together with the analytic 
result in ( |4.1| ). One sees that ( [4.1| ) indeed corresponds to the lowest eigenvalue. 

To see how the approximation by the N c — > 00 result works, we compare the two results 
for n = 4, as an example. The exact mixing matrix is 



Y 



511 

4J 
LO 

3 



45 
10 



21 
122 ) 



(4.6) 



10 



and the one at N r — > oo is 



Y 



119 
" 10 

3 



•14 



(4.7) 



From these Ys' one gets 
M 4 \e(Q 2 ] 



(0.983a 4 , 2 (/u) + 0.520a 4 , 3 (/x)) 



+ (0.017a 4 , 2 (A0 - 0.020a 4 , 3 (//)) 



'a(Q 



2 n\ 9.59/f.o 



'a(Q 



2\\ 15.3/fco 



a(/i 2 



(4- 



for the exact Q 2 evolution and 



Mi 



e(Q 2 



04,20*) + 



2 U 10.4/fio 



(4.9) 



at AT C — > 00. Equations ( |4.8|) and Q4.9Q can be compared with Eqs. (20) and (21) of for 
hi- One sees here again that the coefficients in the second term of ( f4.8| ) are small (i.e., 1/N 2 
suppressed) and the anomalous dimension in (|4.9| ) is close to the smaller one in (|4.8|) . 

5 Summary 

In this paper, we presented the covariant calculation of the anomalous dimension matrix for 
the chiral-odd twist-3 distribution e(x,Q 2 ). The operator mixing with the EOM operator 
is taken into account. This study completes the whole list of the anomalous dimensions of 
all the twist-3 distributions in the one-loop level together with the known results for g 2 and 
h]j. At large N c , we have confirmed that the lowest eigenvalues of the anomalous dimension 
matrix take a simple analytic form shown in Q4.1p and it governs the whole Q 2 evolution of 
e(x, Q 2 ) as was found in [0 by a different method. 
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In this appendix we present the contribution from each one-loop Feynman diagram shown 
in Figs. 2 and 3. 

First we consider the three-point function FyVt^ with the insertion of R U) i ■ A. We set 
m q = to get Zi m and Z iE (Z, m = 2, . . . [n/2] + 1). 
Fig. 2(a) gives 0. 
Fig. 2(b) gives 

,2 



g 



16tt 2 e 



a 



G 



/ (l + m)(m- 1) (n- l + 2 + m)(m- 1) 
^ 2 \2(l-m)[l- 1] 2 + 2(n-Z + 2-m)[n-Z + l] 2/ 



7^(3) 



1 — S;_l — S". 



1 



1 



n-l+l 



+ 



(n + 2)(Z- 1) 



2Z 2{n-l + 2) 2(n-2Z + 2)[n-Z + l] 2y 
+ ( ^ } /(2n-Z-m + 4)(n-ra + l) + (n - Z + 2 + m)(m - 1) 

m = / + l V 



2(m - Z)[n - Z + l] a 2(n - I + 2 - m)[n - Z + 1] ; 



y?(3) 



ft(3) 



1 i n ( 3n -2Z + 6) ^( 3 ) 

+ -t 1 ; t„ » x , n: even 

4[n-Z + l] 2 n-2l + 2 n '2 +1 



2 -''*•'< 5 



(A.l) 



where {■ • • n : even} means this term is only for an even n and its coefficient is one half of 
that of the third line with m — [n/2] + 1 (the same in ( |A.3| ), ( A.5 ) and ( |A.7| ) below), and 

,2 



9 



lQiT 2 e 



a 



G 



E hr 

m=2 \ 2 



n + m + 1 



— m + 1 



7^(3) 



1 - 5" - 5a+i) 



for / = [—] + ! with even n ) . 



(A.2) 



In ( |A.2|) , the terms with m = 2,..., [n/2] are the same as the first line of ( A..1Q with 
/ = [n/2] + 1. This rule also applies to the results for other diagrams. 
Fig.2(c) + Fig.2(d) gives 



lQir 2 e 



(2C F - C G ) 



E 2 (- 1 ) J 



.m=2 



-1 [ n-l+lCm—1 1-lCm-l 

{[n-l + l], J^, 



+2(-l)' 



-1 f n-i+lCn + ( — I)' 



ft 



(3) 



[n-l + l], ' [l-l] 3 J 



{n/2) 2 (-l) m " 1 



(n-i+l^m-l + ( — 1)" n-l+lCm-i) 



2(-l)f (3) 
+ [n-Z + l] 3 ^f +1 " : 6Ven 



1 



1 



[Z] 2 [n-Z + 2] 5 



2</<(| 



(A.3) 
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and 



16n 2 e 



(2C F - C G ) 



[n/2]+l (~< 



+ 1 with even n 



- + 1 

2 ~ x 



(for I = 


~n~ 




.2. 



(A.4) 



Fig.2(e) + Fig.2(f) gives 



16tt% 



(-l)'-™,_ 2 Ci 



i 



l—m 



m=2 



- c G - + 



J n-l+2 
+ {{2C F -C G " 



(Z - m) 

7?(3) 



+ 



n—l—m 



id 



n—l^m—2 



[n — I + 2 — m) 



n— m+1 



J n-l+2 

( (n/2) 

+ 2 (2Cf-C7 ; 

[m=/+l 

-C G ^-—1^11 

ra — Z + 2 



n-2Z + 2)(n-Z + l) 



*2 



2Cp\Si-i + SVi-2+1 — 2) 



-1) 



rn—l r~f 

n—V^m—l 



+ 



i \n—l—m ft 
~ L ) n-l^m-2 



m-Z) m _iC m _j (n - Z + 2 - 77i) n _ m+ iCi_i 



+ (2c F - a 



G 



■1)5- i+1 „-iC 5 _ i+ i' 



and 



Z n-Z + 2 



n/2 



77 — Z + 1 »C». 
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Fig. 2(g) + Fig.2(h) gives 
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Next we calculate Zie, Zin, Zee and Znn from the one- loop two-point functions shown 
in Fig. 3 with a nonzero quark mass. 

The one-loop correction to the two-point function with R n j ■ A comes from Fig. 3 (b) 
and (c). It gives 
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For the one- loop correction to the two-point function with E n ■ A, Fig. 3 (a) gives 
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and Fig.3(b)+Fig.3(c) gives 
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Although each contribution from Fig. 3 (a) and (b) is different from the one for Hl, the sum 
of ( |A.10| ) and ( |A.11| ) are the same as the one for h^. 

Finally, for the one-loop correction to the two-point function with N n ■ A, Fig. 3 (a) gives 



92 -Ce—^UP., 



16n 2 e~ r n(n + l)" n ' 
and each of Figs. 3 (b) and (c) gives the same contribution, 
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These are the same as the ones for hi- 
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(A.12) 



(A.13) 
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Figure captions 



Fig. 1. (a) Three-point basic vertex for Rnj, E, and N. (b) Four-point basic vertex for R n i 
necessary for the calculation of the diagrams shown in Fig. 2. 

Fig. 2. One-particle-irreducible diagrams for the one-loop correction to F^p, q, k). 

Fig. 3. One-loop corrections to the two-point functions. 

Fig. 4. (a) Complete spectrum of the exact eigenvalues of the anomalous dimension matrix for 
e(x, Q 2 ) together with those for hi (squares), (b) Complete spectrum of the leading 7V C 
eigenvalues of the anomalous dimension matrix for e(x,Q 2 ). Solid line is the analytic 
solution in (|4.1| ). 
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(g) (h) 
Fig.2 




Fig.4 
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Fig.3 



